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We show that Inflation in a False Vacuum becomes viable in the presence of a spectator scalar field 
non minimally coupled to gravity. The field is unstable in this background, it grows exponentially 
and slows down the pure de Sitter phase itself, allowing then fast tunneling to a true vacuum. We 
compute the constraint from graceful exit through bubble nucleation and the spectrum of cosmo- 
logical perturbations. 
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I. INTRODUCTION 



Observations in the framework of the standard cosmological model tell us that the observed Universe is in a zero- 
energy vacuum state (except possibly for a small cosmological constant which becomes relevant only at late times). 
" From the point of view of fundamental physics, however, one can wonder why the vacuum energy is zero. The 
^ ' question that we want to answer in this paper is the following: is it possible that the Universe started in a different 
vacuum state with large zero-point energy and then tunneled to the zero-energy state (which could very reasonably 
0^ the minimum available energy state)? * 

' At the same time it is also well-known that a Universe with a large vacuum energy inflates. And, in fact, this 
was the first way to introduce Inflation ("Old Inflation", Guth, 1982 3]). As Guth himself already realized in the 
original paper, the simplest version of this idea does not work, since inflation does not end successfully ("graceful 
. exit" problem). 

A mechanism that could realize this idea (often called False Vacuum Inflation or also First-Order Inflation, since it 
would end through a first order transition with nucleation of bubbles of true vacuum) in a simple way, appears very 
attractive: in fact, it would explain dynamically why the Universe has almost zero cosmological constant, and at the 



O ' same time it would provide Infiation, without the need for a slow-rolling field. In this paper we show that all this is 



possible with the addition of one simple ingredient: the presence of a non minimally coupled scalar field. 



^ This model has already been proposed by Dolgov ^] in 1983 in order to provide a solution to the cosmological 

. . , constant problem: the scalar field is able to screen the bare value of A at late times. However the model was not 
^ ' a realistic description of our Universe, since at late times the Newton constant was driven to unacceptably small 
k>( \ values. Our idea instead is to use this model in order to explain how a False Vacuum Inflation can end successfully 
. through Bubble Nucleation. A model in the very same spirit was proposed already in the '80s, known as Extended 
C3 ■ Inflation (EI) 0,0. The problem with Old Inflation was that, in order to have sufficient inflation one needs the ratio 
r = Tvac/Hf (where Fvac is the decay rate per unit volume of the true vacuum to the false vacuum and Hj is the 
Hubble constant during infiation) to be much smaller than 1. At the same time Inflation can end only if r becomes of 
order one. So only a model with variable (increasing) r can be viable. The idea of Extended Inflation was to consider 
a Brans-Dicke gravitational theory, in which the vacuum energy does not produce de Sitter Inflation but Power-Law 
Inflation. This gives indeed an increasing r and so successful completion of the phase transition. However, the model 
had a definite prediction: the spectral index of cosmological perturbations had to be roughly ns ^ 0.8, and this was 
ruled out in 1992 by the Cosmic Background Explorer satellite (COBE) experiment on the anisotropics the Cosmic 
Microwave Background 7], so many variants of this model were proposed 0,0]. Another possibility for making r 
variable that has been explored is having Fvac variable through the introduction of a slowly rolling field jl^. Our 



* Electronic address: dimarco@bo.infn. it" 
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* It is even more appealing if we follow Coleman and De Luccia arguments 0| , who concluded that in the thin wall limit and for large 
enough bubbles, tunneling is suppressed to a negative energy state, or Banks Q who claims that such a tunneling cannot occur since 
instantons that interpolate to a negative energy vacuum simply do not exist in a consistent theory. 
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proposal still keeps Fvac constant, but it differs from the so-called hyperxetended ^1, P| models, since we introduce an 
initial de Sitter stage and then we show that this stage dynamically evolves in a second stage, which is the same as 
in the original Extended Inflation. In this sense our proposal is as simple as the original EI: the difference is that it 
has two different periods of inflation. 

The first period is exactly dc Sitter and it is the period in which the spectrum of perturbations is produced on the 
scales relevant to the observations (surprisingly none of the modifications of Extended Inflation has considered the 
simple possibility of having an early de Sitter stage, the only mention being in lT|, which however did not consider 
it as a viable opportunity) . During this period the non minimally coupled scalar field grows exponentially from small 
to large values. The second period instead starts when the non minimally coupled scalar becomes important, it slows 
down inflation dramatically and it leads to a power-law behaviour (which could even be decelerating). So it can lead 
to a successful transition to the radiation era through bubble nucleation. 

The plan of the paper is as follows. In the next section we introduce the model proposed by Dolgov Q, summarizing 
the dynamics of the classical homogeneous fields. In section Hill we analyze the constraints on the parameters of the 
model in order to give sufficient inflation and graceful exit, both analytically and numerically. In section ITVl we study 
the production (spectral index and amplitude) of density fluctuations in this minimal scenario, and we also discuss 
the fluctuations due to an additional minimally coupled scalar field which may act as a curvaton. In sectionlVl instead 
we show possible ways to recover Einstein gravity at late times. Section IVII contains a discussion of other possible 
effects and of other observables still to be computed, that are not covered in this paper. Finally, in section IVlII we 
draw our conclusions. The paper contains also two appendices. Appendix ^ contains the calculation of the constraint 
coming from the absence of detection of big bubbles in the CMB. Appendix iBlcontains some calculations of the initial 
amplitude of the field </> and of the perturbation in the energy density due to fluctuations in the non minimally coupled 
field. 



II. THE MODEL 



As we have mentioned in the Introduction, the model that we are going to study has an unstable vacuum energy, 
that we will call A, with some tunneling rate per unit volume Fvac^- 

The only other ingredient of our model at this point is, as in 0, a scalar field cf) which has a generic ("non-minimal") 
coupling to the Ricci scalar R. So, the action is: 



(1) 



where it is crucial to assume /3 > 0. The reader should note that the field is a spectator and it has nothing to do 
with the false vacuum energy A. 

Note also that the bare mass in the Lagrangian M does not coincide in general with the value of the Planck mass 
that we observe today (Mpi). The Einstein equations that follow from this action^ in the homogeneous case are: 
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H ^ -- 



3(M2-h/3</)2) 
1 



2(M2-h/302) 

where H = a/a (a is the scale factor), and the Ricci scalar is: 

and the evolution equation for the scalar field is: 

4> + 3H(p - 13R(I) = 



6i//3(/)0 -I- A 



(2) 
(3) 

(4) 
(5) 



t Note that, although we are using the name A, this is not a cosmological constant, but a vacuum energy that can decay with some rate 
Fvac I as happens for example for a scalar field trapped in false vacuum. This can be in principle computed by knowing the details of 
the potential barrier that separates the two vacua. 

-f In this paper we stick always to the Jordan frame. For the analysis in the Einstein frame see the case n = 2of Ea. 
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FIG. 1; The evolution of the Hubble parameter as a function of time in our model (solid black line), compared to a power-law 
evolution H = a/t (red dashed line). Here /? = 1/56, Hi = IQ~^M. 



The energy density and pressure for the field (j) are: 



1 -^^ 



6H(3(l)(j> - ijSH 



2 J,2 



1 



(6) 
(7) 



Now, let us assume that the Universe starts in the false vacuum, and that the field 4> sits close to zero at the 
beginning. So, at initial time the evolution of the scale factor is given by a de Sitter phase: 

The equation of motion for in such a background is: 

ij) + - l2Hj(3(f) = , 

and its growing solution is: 



e = 3 -1 



(e ~ 8/3 for smaU (3 ) . 



(8) 
(9) 

(10) 



As we said, we are going to choose the initial condition (pQ close to zero (and zero initial 0o). Actually it is not physical 
to assume that (p starts exactly in zero, since the field is subject to quantum fluctuations, and so the initial value is 
of the order (/jq « Hi (see Eq. (jB4ll for a more precise estimate). So the field (j) is unstable in this background, and 
indeed it drives the de Sitter phase itself to an end. In fact, the energy density related to the field (f) is negative: 

p^~~AHUle'"'\ A = 3/3(l + e)-ie2 > 0, (11) 

and increases with the time in magnitude. 

After a period of transition, at late times the field (t){t) becomes linear in t, its energy density becomes constant, 
and the scale-factor increases with time as a power-law: 



where: 



a[t) - , 
_ 1 + 2^ 



B = 



4/3 v/60/32 28/3 3 

A numerical solution of the system is presented in FIC^^- 



(12) 
(13) 



§ Some readers might feel uncomfortable with 4> becoming bigger than the fundamental scale M, being worried about quantum gravitational 
corrections. However these are well known to be under control 13| if the energy density and the masses are below M, which is the case 
here. 
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This late time behaviour is due to the fact that, when becomes very large, the initial "bare" mass M becomes 
subdominant, and so the system reduces to a pure Brans-Dicke model of gravity. It is known that vacuum energy 
in Brans-Dicke theory leads to power-law expansion: this was in fact used as a proposal for the graceful exit in the 
Extended Inflation scenario. The crucial difference here is that we have a dynamical transition between Einstein 
gravity (that leads to the exponential phase of inflation) and the Brans-Dicke gravity (that leads to the power law 
phase of inflation) . We will exploit the presence of the first phase in order to produce a flat spectrum of cosmological 
perturbations, and the presence of the second phase in order to have a graceful exit from inflation to the radiation era: 
in fact H decreases and when it becomes of the order of the decay rate r^ac it allows the system to tunnel successfully 
to the true vacuum. 



III. CONSTRAINTS ON INFLATION 



We show in this section under which conditions the model provides a "graceful exit" from inflation. The transition 
from the "false" vacuum to the "true" vacuum with zero energy density happens through bubble nucleation. The 
transition has to be sufficiently abrupt, so that essentially all of the bubbles are created in a very short time. The 
nucleation in the early stages of inflation, in fact, would produce bubbles that would be stretched to very large scales 
by the subsequent inflationary phase, and that would spoil the observed isotropy of the CMB on large scales [l^lTELIl^ 
("Big Bubbles" constraint). 

First, we introduce the following definition for the number of e-folds: 

N{t) ^ J^""" H{i)di, (14) 

where tend is the time at which inflation ends through bubble nucleation and the radiation era starts. Assuming that 
inflation starts at t = 0, we define A'tot = ^(0) as the total number of e-folds in inflationary stage. Our model can 
be successful if it is able to produce the correct spectrum of perturbations on large scales and, as it is well-known, 
a flat spectrum can be achieved if H is sufficiently constant during inflation. This seems to be in contrast with the 
requirement that in our model H has to decrease in order to allow the phase transition. However, the crucial point 
that we want to stress here is that the requirement that we need in order to satisfy observations is not so strong: we 
only know that the scales between 3000ft-~^Mpc and about 50ft.~^Mpc have a flat spectrum This can be achieved 
if these scales go out of the horizon during the first stage of inflation, in which the Hubble parameter is constant. In 
usual inflationary scenarios, these scales are produced roughly between 60 ^ ^ 64. For the following we define 
Nl = N{tL), where is the time at which a given scale L crosses outside the horizon. We also define the phase 
I as the phase in which H ~ Hi is a constant, going from the initial time t = up the time ti where a transition 
happens. In terms of e-folds it goes on for iV/ = iVtot ^ N{tj) . The phase II starts at tn (slightly after i/), when H 
begins to evolve such a/t, and ends when the vacuum decays, at a time tend, and goes on for Nn = N{tij) e-folds. 
There are three basic requirements that we must satisfy: 



• on the relevant scales {i.e. during the phase I ) the production of bubbles of true vacuum has to be sup- 
pressed, and the only source of perturbations have to be the usual quantum fluctuations of scalar fields, 

• the phase II has to start after the production of perturbations on these scales: Njj ^ AsQ/j-iMpci 

• inflation has to be sufficiently long: TVtot ?i A^3ooo/i-iMpc- 

The crucial quantity that regulates the production of bubbles is: 

^<"-W?^ (15) 

The end of inflation tend is achieved when this quantity is of order 1. More precisely f^e condition for the 
percolation of bubbles is: 

r(tend) = ^ ■ (16) 

We show in FIG. |2] the typical behaviour of r{N) which is clearly constant for large N and it increases with 
decreasing up to ^ (where A^ is defined to be zero). 

The value of the ratio during the phase / is a free parameter, but it is constrained to be small: 

Fv 



ro = ^<10-^ (17) 



Hf 
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FIG. 2: The evolution of r = V^^^/H'' as a function of the number of e-folds N (with 13 = 1/56, Hi = 10"^M, ro = 10"'^). 
Time goes from right to left: r = ro at the beginning and then it grows to 0(1). 



in order to prevent big bubbles formation. As we mentioned, in fact, a bubble nucleated so early would have very 
big size after inflation and it would appear as a black patch in the microwave sky. The number given as a bound is 
computed in Appendix 

So, what is important is to constrain the ratio r[N) as a function of the number of e-folds. 

In order to give a quantitative constraint we need the correspondence between a scale L and the N at which they 
are produced (A^l), which is given in general by: 



L 



L ' 



(18) 



where is the value of the Hubble constant when the scale L exits the horizon, and Trh is the temperature at the 
beginning of the radiation era. So, the relation is: 



Nl = 63.3 + AiV + In 



3000 h-mpc 



(19) 



Since we assume that the scales between 3000/i ^Mpc and about 50h ^Mpc are produced during the phase I, we get 
the following two relations: 



iV30oo?>-iMpc - (63.3 + AiV) , Wsoh-iMpc - (59 + AiV) 



(20) 



where: 



AN = 3.7 + In 



Hi 



(21) 



In our model we may assume that reheating is almost istantaneous due to rapid bubble collisions, and so: 



1/2 r72/,^(R)N2 



(22) 



where the superscript "R" means that it is the effective value of the Planck mass at the beginning of the radiation 
era. This implies: 



AN ~-\n 
2 



H, 



lO^GeV 



1 



In 



10- 



1 



■In 



(R) 
PI 



lOi^GeV 



(23) 



Since the amplitude of the cosmological perturbations constrains Hi to be at most 10 (see Eq. ()55|l). the biggest 
value that we can get for AA^ is around: 



AN 



-lO-' 



(24) 
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A. Analytical approximation 

We give here an understanding of all the constraints showing their dependence on the parameters of the model with 
an analytical approximation. A more precise numerical analysis of the ratio r{N) is provided numerically in the next 
subsection. The number of e-folds of the first phase is given by: 



Ni ~ Hiti 



(25) 



where the time ti of the end of the first phase is found equating the energy density in the field (p to the bare vacuum 
energy A: 



A ~ \p^{ti)\ = AHj^itiY 



Using Eq. (|10() we have: 

and so the number of efolds is: 

Ni 



1 



V9T48/3- 3 
It is of some utility to expand it for small (3: 



■In 



A 



02 (V9 + 48/? (1 + 4/3) - 3 - 16/3) 



iV..3^(21n 



M 
4>a 



- ln[/3] 



(26) 



(27) 



(28) 



(29) 



It is in principle more difficult to estimate the duration of the remaining part of inflation, since we do not have an 
analytical expression for the transition between the phase I and phase II , when H drops from the value Hj to a/t. 
However, starting from the end of inflation we can use the asymptotic solution H = a/t and so we get from Eq. H14|l : 



Nit) = a In 



We have also that in the asymptotic phase the ratio r(N) goes as: 



r{N) 



4 ''end 



a"* 47r 
We can roughly use this solution to find when the phase //begins: 

9 



47r 



that leads to: 



--In 



An 



So, imposing Njj < Nc^^i^- iy[p^ we have a constraint on a (at given tq) 



< 



-4(59 + AA^)/ In 



Now inserting here the maximal value for A A of Eq. H24I) , we obtain: 

-236+lln[^] 
ln[^] 



(30) 



(31) 



(32) 



(33) 



(34) 



(35) 
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This translates in a constraint on f3 (or equivalently on iv = which is sometimes used in the hterature) . Assuming 
the maximal possible value tq ~ 10"^ we get: 

/3>-^ or uj<U.5. (36) 
oo 

As we said this number is only an estimate, and we have the correct constraint coming from the numerical analysis, 
in Eq. (|37|l . It is relevant to stress that this condition is independent on the third parameter of our model, that is 
the value of Hj/M (or equivalently the value of the vacuum energy A). As we will show in section HVI the range of f3 
that we get will severely constrain the spectral index of the fluctuations in the field (j). 

Then, the third constraint that we have to impose is that the duration of inflation is sufficiently long. In terms 
of e-folds, the duration of the first phase is given by Eq. 128(1 . while the duration of the second phase is given by 
Eq. ((33|l . so one can check for a given set of the parameters whether Nj + Nu > -/Vsooo/i-iMpc^ is satisfied or not. 
There are three relevant parameters for this condition: 4>q/M (which is actually Hj/M if the initial condition is 
given by quantum fluctuations), P and tq. The condition is very easy to satisfy for small P and more difficult for f3 
approaching the value 1/2 (for /? > 1/2 we do not have inflation anymore in the phase 11). 

B. Numerical analysis 

In this section we present the result from the numerical solution of the differential equations that govern our system. 
In particular we are interested in r{N) and in constraining it. 

The way in which we applied the conditions on the duration of inflation given above is as follows: 

• we compute r{N) — rvac/^^(A^)* and impose that ^(Afso/i-iMpc) < 10^'', 

• we impose that the variation of H between -/Vso/i-iMpc and -/V3ooo;i-iMpc is small **, 

• we check that iVtot ^ A^3000/i-iMpc- 

The region of parameters in which the model works is shown in FIG. 13 for some different values of AA^. 

As we showed in the previous subsection, there are two main parameters: the ratio Fvac/^^/ and the nonminimal 
coupling parameter /3. As we have shown qualitatively in Eq. I|34|) . the viable interval in the parameter (3 shrinks as 
To decreases. This is respected by the behaviour of FIG.|3| except for the upper part of the figure, where the relevant 
scales are produced exactly at the transition between exponential inflation and power-law, and so r(N) decreases 
slower than our analytical approximation ta (N) of Eq. H31|) . 

The most important information that appears from the exclusion plot is that the mechanism works with an upper 
bound on (3. Considering that AA^ < 0, from the plot we can see that: 

13 >^,. (37) 

which means that the analytical approximation Eq. (|36|l already gave a good estimate. 

Also, note that the dependence on Hj comes only in two places. First, Hj enters in defining what is AA''. Second, 
when checking if the model is able to give sufflcient inflation (which is not shown in FIG. IS)). In fact, the duration 
of the exponential phase is related to the ratio (j)o/M (see Eq. ((28|l '). and (po ~ Hj if we assume that (j>o is set by the 
quantum fluctuations (see Eq. l|B4p ). It becomes difficult to have sufficient inflation only for relatively big values of 
(3 (close to 0(0. 1)) and when Hj is close to the scale M. 

IV. SPECTRUM OF PERTURBATIONS 

In this section we compute the spectrum of cosmological perturbations produced in our scenario, in order to see if 
it is consistent with the observations. 



^ Here we are disregarding the transition between phase I and //, so that A^tot ~ + Nij. 

** By "small" we mean here an amount that does not change substantially the spectral index. Precisely we take 2AH/AN{1/H) < 0.1, 
where AH = ^^CA^aoooh-iMpc) "^(-'^SOh^iMpc)- ^^'^^ ^^'^ ^^'^ correction that would get a spectral index in slow-roll approximation, 
and we want it to be smaller than 0.1. 
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FIG. 3; In grey we show the regions of allowed parameter space, for different AA^ (from left to right AA'^ = —8, —4, 0). The 
region on the top and the region on the right are excluded because of the constraint on the production of large bubbles (this is 
independent on the value of Hi /M). For /3 becoming bigger than what is shown in the plot, the scenario is still viable except 
when P gets closer to 0(0.1), where sufficient inflation might become difficult to obtain (depending also on the value of Hi /M). 

The perturbations on relevant cosmological scales are produced in our model during the phase of exponential 
expansion. For this reason one can realize situations in which the observed spectrum is flat, i.e. ns = 1- 

For comparison Extended Inflation instead has always a power law expansion , where the Big Bubbles 
constraint requires a < 20 as we discussed in sect. IIIII [T^. The constraint on a produced a spectral index ns ^ 0.75 
for any fluctuating field, which was clearly ruled out by observations, after COBE 0- Even adding a curvaton did 
not help p^ . since any field has generically a non-flat spectrum, due to the fact that H is never constant. 

In our model, first we analyze the slope and the amplitude of the spectrum of perturbations due to fluctuations of 
the field cj). Physically the fluctuations in (p (which is responsible for ending inflation) give rise to end of inflation at 
different times in different Hubble patches, therefore providing a mechanism for density perturbations. However, even 
with the presence of our exponential phase, this minimal scenario appears to be in disagreement with the observed 
spectrum. The reason is that the field (j) is not minimally coupled and so its fluctuations do not have ns = 1 even if 
the Hubble constant is exactly constant during the first phase of inflation. As we will show the situation is improved 
with respect to Extended Inflation (precisely — 1 is a factor of 2 closer to 1), but the Big Bubbles constraint still 
constrains ns to be too small, using the WMAP data. 

One possibility (which is under study and it is not developed in this paper, except for some remarks in sect. lVI|l 
to construct an inflationary model fully compatible with observations is to modify the original Lagrangian in order 
to obtain a strong successful transition even with a small /3 during the phase I. Another possibility is to keep our 
Lagrangian, noting that a good spectrum can easily be produced assuming that another light minimally coupled field 
(so, with flat spectrum of fluctuations during the exponential phase) is responsible for the generation fo the curvature 
perturbation. We will show the viability of this possibility in subsection IIVBI 



In order to get the spectrum of the cosmological perturbations we consider fluctuations 5(j){x, t) associated to the 
scalar field </>: 



and we follow the treatment of cosmological perturbations in generalized gravity theories given by |20|, |2l| . In a flat 
cosmological model with only scalar- type perturbations: 



A. Computing spectral index and amplitude with only (j) 



(t>{x,t) = (j){t) +6(l){x,t) , 



(38) 



ds^ = -(i + 2^/j)dt'^ - x,^ dtdx^ + a2(l + 2^)5,,^ dx^dx' 



(39) 
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we may study the evolution of the fluctuations by means of the gauge-invariant combination: 

6<j)^ = 5(1)- j^ip, 

which is 5(f) in the uniform-curvature gauge ((^ = 0). The equation of motion for 5(j)^ is: 



(40) 



H 

a^Z^ 



a^Z 



H 



where the function Z is defined as: 



Z{t) = 



+ /302)(M2 + /3(/)2 + 6/32 ,/,2) 

{AP + f3(P^+f3(t>{^/H))^ 



In the large-scale limit the solution is: 



H 



C{x) - D{x) 



a^Z(j)^ 



(41) 



(42) 



(43) 



and, ignoring the decaying mode D{x), we can relate the spectrum of 5(j)^ to the time independent function C{x) 
which represents the perturbed 3-space curvature in the uniform curvature gauge. So, we can compute the amplitude 
of scalar perturbations by means of the amplitude of the variable C{x). 

The next step is to rewrite Eq. (|41|l (in fc-space and in conformal time r) in terms of the variable v = a^/Z5(j)^, and 
we can take into account that, during the phase I (in which the perturbations are produced), the function Z reduces 
to Z ~ 1. So the equation of motion becomes: 



(44) 



This result agrees with the general fact that a conformally coupled field (/? = — i) has only Minkowski space fluctua- 
tions. We can define the conventional constant v as: 



which, for small /3, becomes: 



The solutions of Eq. iBUl are: 



v = - + - , 
2 2' 



v{k,r) = -^JAA [ci(fc)i/(i)(fc, |t|) + C2{k)Hl^\k, \t\) 



(45) 



(46) 



(47) 



where ci(fc) and C2(fc) are constrained by the quantization conditions |c2(fc)|2 — |ci(fc)|2 = 1 and ff^'^ are the Hankel 
functions of the first and second kind. So, on large scales the power-spectrum of C is: 



H 



Psilik,r) = 



H 


H Til.) 


■f k ' 


J 


2-r(i) 


2^ 



|ci(fc)-C2(A:)| 



(48) 



^From here we can immediately extract the spectral index ns for scalar perturbations (using the Bunch-Davies vacuum 
C2 = 1): 



ns = l + -— ^ - 4 - 2j/ = 1 - e , 
dink 



(49) 



and for small /3 this becomes: 



ns - 1 - 8/3 . 



(50) 
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As we said, we have /3 > 1/57, so ns ^ 0.86. Unfortunately this is not in agreement with recent observations (such as 
Wilkinson Microwave Anisotropy Probe (WMAP) [13]); therefore the model as we presented requires some corrections 
in order to work, as we said, modifying the Lagrangian or adding a curvaton. 

Let us compute the amplitude of the perturbations in 5(j)^ so that we will impose that it is suppressed (in case we 
add a curvaton this is required). The amplitude of the spectrum is: 



Ac 



Hi T{u) 
^e(/.r(|) 



(51) 



This has to be evaluated at A'^3qqq^-ij^p(., therefore (f) at that time has to be much larger than _ff/, if we want the 
amplitude of the perturbations to be small, as required by experiments. Since we are assuming the initial value (/jq 
to be of order Hi , this means that there must be a sufficient number of efolds between the inital time and the time 
corresponding to A^3ooo/i-iMpc- This can be achieved if the value of 4> at the end of the phase I (which is roughly 
M/v^) is much larger than Hi. Therefore this tells us that M has to be large enough with respect to Hi. Precisely 
we may write (f){t) ~ M / y^e'^^^^*'^''' (valid for small /?) and so the amplitude for the mode corresponding to present 
horizon, reduces to: 

Ac = "I .8/3Ajv,„, (52) 
64/37r2M2 ^ ' 

where AiVhor = A^soooft-i Mpc - N{ti). 

In order to connect to constraints form observations let us also introduce the total perturbation in the energy 
density C, which is gauge-invariant and defined as: 

(53) 

Ptot 

where ptot is the total energy density, and (5/Otot stands for a perturbation about the average value. One can compute 
it in the uniform gauge and find that ({x) = C{x) neglecting the decaying mode in the Eq. (|43|l (see appendix 01. 
The quantity (^{x) is conserved on superhorizon scales in any cosmological epoch, and therefore it is related to density 
perturbations at late times and it is constrained by experiments. Precisely, its amplitude has to be smaller than 
4 X 10~^°, if we want to suppress the spectrum (since its spectral index is not flat). So this requirement constrains 
the ratio Hi /M (we give here the result for small (3) : 

^ < le^v^e-^^^^-' X 10-5 . (54) 

This condition guarantees a sufficient number of e-folds from the initial time Iq until the time at which the scale 
corresponding to the present horizon is produced. 

We can also work out a necessary condition: since the right hand side has a maximum for {3 = l/(8AAfhor) and 
ATVhor > 4.3 we find: 

^<5xl0-^ (55) 
which is also necessary to satisfy the phenomenological bounds on tensor perturbations (see Eq. Hfi5(l ). 



B. Adding a curvaton 



As we said at the beginning of this section, another field (cr) may be responsible of the generation of the curvature 
perturbations. This is necessary (unless we modify the Lagrangian to have a strong transition even when /3 is small 
during the phase I, as in ) since the spectral index of perturbations in (j) is not close enough to 1 . In order to do 
that we require the field a to be minimally coupled (or with very small nonminimal coupling) and light. First of all, 
note that if Hi happens to be small enough (see Eq. 1)54(1 ^ then the amplitude of the fluctuations generated by the 
field (j) alone is negligible. In this case it is the field a that can give rise to the observed spectrum of perturbations. 



tt It is immediate from the second of our basic requirements in sect. If If I and from Eqs. l2Ui . 
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through the curvaton mechanism |22| . So, we consider here a minimally coupled field cr, whose mass m„ is smaller 
than Hj. In other words we have to add to the initial lagrangian of Eq. the following terms: 

Lm^-^d^ad'^a-^mla^ (56) 

where for simplicity we have chosen a quadratic potential. Since the field is light it develops quantum fluctuations 
with an almost scale invariant spectrum during the phase I, in which H is constant. The spectral index is given by: 

= (57) 

The reason why it differs from the (j) fluctuations is that fluctuations in a minimally coupled field are sensitive only to 
the slow- roll parameters (that is the evolution of the background) , while a non minimally coupled field is also sensitive 
to the couphng 

Also, since it has no coupling with 0, its energy-momentum tensor is conserved and so the variable: 

C. ^ -i/^ (58) 

Per 

is constant on large scales. In fact, as it was shown in [2^, the quantities: 

Q^-^-H^-£l (59) 
Pi 

(where the subscript "i" refers to one single component of the energy content of the Universe) are conserved on 
large scales if the energy-momentum tensor of the i component is conserved, irrespectively of the gravitational field 
equations. So this applies to generalized metric theories of gravity and so also to our case. Now, the total C can be 
expressed as a combination of the individual Q: 

J. Crest Prcst ~t- C^crPcr /r'r\\ 

C = ■■ —■ , (60) 

Prcst + Pa 

where the subscript "rest" stands for any other component in the Universe. So, if evolves differently from prost, 
then C is no more constant. In particular, if the numerator is dominated by C^Ptn then the variable C acquires also a 
flat spectrum. 

The curvaton scenario uses the fact that when H becomes smaller than m^, the field starts oscillating. In our case 
we also need that ma- ^ tH^, so that the oscillations start after the end of inflation t^;, so everything proceeds as in 
the conventional curvaton scenario. Therefore the "rest" is the usual radiation component during the radiation era. 

In fact, a is frozen to its initial value (cr/) before the time tend (which is the moment at which radiation is created 
by the nucleation and collisions of bubbles) and then it starts oscillating when H ~ rria- So at this point it starts 
behaving like matter and eventually it becomes dominant as in the usual curvaton models |22j |. since it redshifts less 
fast than radiation. Therefore its Co- Pa — ~'iHC,aPa can dominate and so C — Co-- this case the amplitude of the 
power spectrum is given by: 

^1/2 Hi 

V/ =rc , (61) 

^ TTai 

where rc is the ratio of energy density in the curvaton field to the energy density stored in radiation at the epoch 
of curvaton decay (it is constrained to be 1 > rc ^ 10~^), and cr/ is the value of the field a during the phase I of 
inflation. The amphtude of Eq. (|()T|l has to be the observed 10~^ 



tt Another way to see it is to compute the perturbations in the Einstein frame |12| . In this frame the field </< becomes slowly rolling, with 
the e parameter very suppressed in the exponential phase and the -q parameter proportional to /3. So, the spectral index n — 1 for </> is 
proportional to r] and therefore to /3, while for a curvaton field it is suppressed, since it is only proportional to e and not to r) |22| . For a 
more general analysis of perturbations for a two- field system in the Einstein frame, the case of the present lagrangian is also recovered 
by generalized theories of gravity examinated in l23l . 

Note that it is reasonable that C,^ (see Eq. 1511 ') is suppressed with respect to C,^. In fact, since <f>i is unstable, it is natural to have 
<t>l > cr/. 
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In order for this mechanism to work, the curvaton has to dominate the universe before the nucleosynthesis epoch, 
so: 

2 " ' \m„MpJ ~ ^ ' 

Now, the maximal value for is Fvac and the value for ct/ is around 10^77/, so there is a possibility of making the 
mechanism work if: 

iJ/ > 10^GeVro^/^°, (63) 
which generalizes the usual result that Hj ^ lO^GeV in the simplest curvaton scenario [25l |. 

C. Gravitational waves 

Here we stress that in our scenario gravity waves are created in the usual way on cosmological scales during the 
phase of exponential expansion. As long as (f> is small so that we are in the phase I, the spectrum for the gravity 
waves is exactly the usual one and exactly scale invariant: 

This is constrained to be smaller than « 10^ by CMB experiments so: 

< 10"^ , (65) 



M 



which is consistent with Eq. H55|l . 



LATE TIME GRAVITY 



Since the parameter f3 has to be big enough (see Eq. \67}i ). it is apparent that the model does not lead directly to 
Einstein gravity in the late Universe, but to Brans-Dicke gravity. In fact an additional scalar may mediate a "fifth 
force" between bodies at the level of the solar system, spoiling the successful predictions of General Relativity. 

We may distinguish two possibilities. If the scalar field is very massive (with respect to the inverse of the solar 
system lenght) then its influence is negligible in solar system experiments, since it mediates a force with too short 
range. 

If, instead, the scalar mass is small with respect to the inverse solar-system distance, then it must be presently very 
weakly coupled to matter for the model to be consistent with observational data. Deviations from General Relativity 
in scalar-tensor theories are usually parametrized, defining 

F((/.) = (A/2 + /302)^ (66) 

by the following Post-Newtonian parameters: 

. _ {dF/dcj,f _ 1 F{dF/d^) d7 

^ F + 2{dF/d(j)f ' ^ A2F + 3idF/d(j))^ d(j)' ^ ' 



The limits coming form the experimental bounds on these parameters |26) lead to the following constraint: 

< 10"^ . (68) 



F + /302 

This was recognized as a potential problem also in the Extended Inflation scenario, since in the original Lagrangian 
there is no mass for the (j) field. However the problem is not so hard to tackle, since there is a long time evolution 
in the system between the inflationary epoch and the epoch in which gravity is tested. So, just as for the Extended 
Inflationary scenario we can use different possibilities for recovering Einstein gravity at late times. 
There are several strategies to overcome this problem: 



13 



• Drive 4> at late times to a value smaller than M, so that Eq. (|68|l is satisfied for any (3. This is possible for 
example adding a potential term for (f> that drives it to zero. In this case M really corresponds to the observed 
Mp\. In order not to change the previous discussion, the potential V has to satisfy V A. Note that in this 
case we have Einstein gravity irrespectively of the value of /3, and irrespectively of the fact that the field has a 
mass today. 

• In a similar way one can imagine a potential V{(j)) (that again has to satisfy the condition V ^ A ) that locks 
the field (p at some generic value (which can be also bigger than M) giving it a mass bigger than the inverse 
of the solar system lenght. In this case there may well be significant post-inflationary evolution, and the value 
Mpi today may be significantly different from the value M. The fifth force experiment constraints are avoided 
since the field is massive today. 

• Another strategy is to modify the model (as in hyperextended inflation "h], or in l27l|'). in such a way that /3 
is not a constant, but it can vary with time. In this case it is sufficient to have a dynamics such that (3 can end 
up having a very small value at late times, so that Eq. Ht)8|) is satisfied. 

• Finally another interesting strategy is to couple (j) differently to the matter sector (as in ) and to the vacuum 
energy. This is a clear procedure in the so-called Einstein frame (in which gravity is described by the usual 
RMpi term only). In our frame the modification consists in substituting the A term with: 

Sa = -^J d'xV^gA^P^^j . (69) 

The addition of the new parameter 77 (introduced by the generalized Brans-Dicke model of makes possible 
to satisfy all the constraints derived in section ITTll without contradicting current experiments on gravity. 

Models with a potential are simple possibilities, but on the other hand one may argue that they introduce small 
parameters. It is true in fact that, in the sense of 30], the potentials have to be very flat: if one computes the ratio 
AV/{A(p)'^ (where AV is the change in potential energy of the field during the whole duration of inflation, and Aip is 
the variation of the field) one finds that very small values are required. However we think that this is not fine-tuning, 
since the only physical requirement is that V A (where a small hierarchy is sufficient): then the fact that A^** can 
become much bigger than 1^ is a result of the dynamics and it is not put by hand. 



VI. OTHER CONSEQUENCES AND POSSIBLE EFFECTS 

In this section we mention some of the observational consequences of the model that have to be explored and that 
are not covered in this paper. 

One is the production of bubbles on the small scales, which could be a striking signature of the model. If bubbles 
are produced at small scales, these could be detected as additional power in the CMB spectrum at high multipole 
{I ^ 1000), and as presence of many voids on the small scale of the galaxy distributions. Already a few groups have 
investigated in this direction ^31]- 33]. and we will explore in future work the possible imprints of our specific model 
in the observations on small scales. 

A second consequence is that our model can automatically incorporate a change of the spectral index at some small 
scale, which might be necessary to fit the data as suggested by the WMAP In fact the spectral index of any 

fluctuating field becomes smaller at the scale corresponding to the transition form the exponential inflation to the 
power-law inflation. 

Another effect to be computed is the non-gaussianity in the cosmological perturbations produced in this specific 
model: for example in case these are produced by a curvaton this is likely to be relevant. Finally also the production 
of gravity waves from Bubble collisions might be potentially observable. 

Also, we mention here some possible variations on the model and some points that are still missing in the analysis 
of this paper. 

First of all, in the presence of a time dependent background metric, the decay rate Fvac could acquire a time 
dependence |34| . This could lead to a different constraint on the parameters of the theory. In particular the constraint 
on (3, and thus on the spectral index of the (/> fiuctuations can be different (for example an order of magnitude difference 
in Fvac at the beginning of inflation and Fvac during the asymptotic stage leads to a 0.02 difference in the spectral 
index). This effect will be subject to further study. 
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Then, there are ways to make the model viable even with the field 4> alone. One might in fact invoke a different 
coupling of the field <j) with the Ricci scalar. Instead of the (3R4P' coupling one may consider for example: 



K R In 



R 

M2 



(70) 



This was proposed in |35| , where the author suggested that such a term in the Lagrangian can arise through quantum 
corrections. In this case the cosmological dynamics is described by an effective time dependent quantity: 



/3cff = K In 



R 
M2 



(71) 



where /t is a number. This is justified as long as the variation in /?eff is sufficiently slow that its time derivatives may 
be neglected. As R decreases, /3eff incresases thus making the process more efficient. In this way the bound on the 
spectral index of the perturbations of the field (j) gets relaxed. 

The same thing applies for any quantum correction to the parameter /3, that could change the bounds in a relevant 
way. 

Generally speaking our idea of having exponential infiation slowed down by an unstable scalar field could be 
easily implemented in other variants of the model, which make the transition stronger so that the curvaton becomes 
unnecessary Basically wath is needed is to have a coupling with R that grows with 0, faster than the (t)^R 

coupling analyzed in the present paper. 



VII. CONCLUSIONS 



We have shown in this paper that a false vacuum can successfully decay to a true vacuum, producing infiation, in 
the presence of a non minimally coupled scalar </>, since Exponential Infiation is slowed down to Power-Law Infiation. 

Then we have analyzed the constraint coming from the fact that we do not want to produce large bubbles that 
would spoil the CMB. This implies the following constraints: tq = TvacZ-ff/ ^ 10~^ (where Fvac is the tunneling rate 
per unit volume of the false vacuum to the true vacuum, and Hi is the Hubble constant during de Sitter Inflation) 
and the nonminimal coupling has to be /? ^ 1/57. 

If (j) is responsible also for primordial density fluctuations then the spectral index is 715 ~ 1 — 8/3 ^ 0.86, which is 
in disagreement with observations. 

However, if there is also a minimally coupled scalar, it can produce a flat spectrum through the curvaton mechanism. 

Moreover, generally speaking the same idea can easily be adapted modifying slightly the model in order to have a 
stronger transition to Power-Law infiation, and so without the need for a curvaton field. One example [l^ is given 
by generalizing the coupling (jy^R to a generic function f(<p)R, where f{(f>) grows faster than (jp for large (p. Another 
example is to couple a slightly different function of R to (j)'^ (see section IVH) . 

Finally the model can lead to other observable consequences (as discussed in section rvT|) . as a change in the spectral 
index ( "running" spectral index) at some small scale , the production of bubbles on the small scales (detectable as 
voids in the large scale structure and at large / in the CMB) and possibly the presence of some non-gaussianity. All 
these effects are subject to future work. 
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APPENDIX A: BUBBLES ON THE CMB 



Here we follow to estimate the constraint on — T^ac/Hj, using the WMAP data. 

As we have mentioned, the presence of a bubble would be directely detected as a void region (so a region with 
Sp/p = — 1) in the CMB, if the bubble size (whose comoving value will be called r) is large enough with respect to 
the resolution (corresponding to a comoving size rf,) of the experiment. The fluctuation in the temperature (through 
a fluctuation in the Newtonian gravitational potential <&) that it would produce is: 



ST $ 1 fSp 



T 3 3 \ p / \L\ss 



r 



(Al) 
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where 5p/ p = —1 for voids, and where Li^s is the comoving horizon size at last scattering (and the Doppler contribution 
has been neglected as in [ijj)- If the bubble size is smaller than the resolution, it can still be detected, but the 6T/T 
is reduced by a factor {r/rbY'. Also, if a bubble is smaller than the last scattering surface width (whose comoving 
value we call r^), it is reduced by a factor {r/r^)'^- So, the minimal bubble size that can be detected is: 

« (3 X IQ-^Ltrlrl) « 6Mpc , (A2) 

where we have used the following values: L\ss — 190ft.^^Mpc, rt, « 7/i^^Mpc and — 6.7/i~^Mpc. Any bubble which 
produces a fluctuation bigger than what is observed (some 10"^) is in contradiction with observations. 

Note that the pixel size of the WMAP experiment is 30 times smaller than the COBE's 17], but we have used only 
the data up to Z w 600 (which corresponds to a comoving scale of about 7/i^^Mpc), for which the signal-to-noise ratio 
is less than 1. Such a small resolution is what makes the bound stronger today, with respect to the numbers used 
byH. ^ 

The requirement in order to be safe from seeing big bubbles is that we demand that the number of voids inside the 
horizon be less than that which gives a 95% confidence level that at least one is in the last scattering surface. This is 
true if: 



1 AttL^t /dNf 



3 y^,^ {At:/3)L% V dr 



dr <1, (A3) 



where ^^^f- is the spectrum of bubbles {Nb is the number of bubbles) which is generated by a specific model and 



dr 



where Lh — 3000/i^ Mpc is the present horizon distance. 
This spectrum, in our model, is calculated as follows. 

A bubble nucleated at time ti with zero initial size grows at the speed of light within the expanding universe and 
it has size x at time t: 



* dt 



x{tM)^a{t) / — . (A4) 



At time ti the bubble nucleation probability per unit comoving volume is 

dNB 



=rvaca(ii)'. (A5) 
dti 

We can change variables to x and evaluate the spectrum at the end of inflation (iond): 

dNB _ ^3 

^-r_a(te„d) (A6) 
and then, compute it in a sphere whose radius Ih corresponds to the physical size of our present horizon at that time 

dNB rr4 47r 3 ajh)^ dh 

-^-r^Hj—lH^^^—^^. (A7) 

The evolution of the scale factor a{t) (and so the ratio ^) could be evaluated numerically. However, we can use a 
rough approximation that already gives us a correct estimate: we assume that H[t) is constant until it reaches the 
value a/t and then the scale factor goes as a power law i". Using this we obtain the following spectrum: 

dNB ^ 47rro HjHl^lj,{a - l)^ 

dx 3 ((iJ//i7end)"-i?/(a-l)ifc„dx))4 ' ^ ^ 

where we considered values of x that are much bigger than the horizon at the end of inflation (i/end)- 

Then we impose a correction factor {g) due to the fact that the voids expand faster than the background, and in 
principle there is also a correction due to the fact that voids get filled by relativistic matter (so their size is reduced 
by a lenght A), so the true comoving lenght of a bubble is: 

r = 5 (^x - a) . (A9) 



H 
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However, in a Cold Dark Matter scenario matter becomes nonrelativistic very early, so the void filling lenght A is 
negligible. On the other hand, in a CDM dominated universe the void growth is non-negligible pjj and it is: 

g ~ 1.85/1^/^ . (AlO) 
At this point we can extract the quantity in which we are interested: 



(INb 47rro 1 / Lnia-l) 



Note now that the quantity IrHe is very big: 



(All) 



so we may safely neglect the first term in the denominator and we get the very simple result: 

cLNb _ 47rro g^L]j 



(A13) 



dr 3 r4 

Finally, by integrating this in Eq. (|A3|) we obtain the following: 

ro < 10"^ . (A14) 

APPENDIX B: AMPLITUDE OF 00 AND C 

It is interesting to understand on physical grounds what should be the value for 00- If we put for example exactly 
00 = as initial value (so that classically = would be an (unstable) solution of the equations of motions at all 
times), we would find that the amplitude of the fluctuations of C diverges for 0o — 0. However the fluctuation 5(j) is 
not divergent, and so the value that 5(j) assumes in one Hubble patch should be included into the initial classical value 
(500 • In other words the minimal value of 0o should be taken as given by the typical value of the quantum fluctuation 
of the fleld. This value can be taken from the expression: 

As fcinin we choose fcjnin = oo-ff, the wavelenght produced at the beginning of the inflationary phase. As fcmax we 
choose instead fcmin = a(ti)H, the last produced wavelenght at the time deflned by a(ti). As a result, we find: 



\Scf>{x)\'^^\r{3/2 + e/2)f — 



- 1 



(B2) 



Substituting ^ = g^"""^' we get: 



\S^{x)\ = ■^r{3/2 + e/2)-^\/e(^'--^')^-l. (B3) 

This means that the field stays constant for about 2/e e- folds and then it starts to grow. So the physical initial value 
that the field takes (assuming that classically the field is initially set to zero) is the value of the quantum fluctuations 
after about 2/e e-folds from the beginning. This means that the minimal initial condition is given by: 

^^^2Hj^^Wl^lI^ (^^^^^Hj forsmall/?). (B4) 

Going back to the computation of the amplitude of the perturbations, the last step consist in relating the spectrum 
of C to a conserved quantity during the post inflationary evolution. One possibility is to use the variable deflned in 
Eq. (58) of which is conserved superhorizon and it is equal to C. 
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Another possibility is to use the variable C that we have defined in Eq. (|53|l . where we may insert Sptot and ptot 
following ,2^] . In particular (5ptot corresponds to Eq. (33) of where one has to insert the solutions of the system 
of Eqs. (54-56) of that work. In fact the energy density and pressure defined in Eqs. (30, 31) coincide with ptot and 
Ptot and they are conserved, if one sets the potential V defined in those equations equal to a constant A: 

(B5) 
(B6) 

where F = M"^ + f^cj)^ . One can verify explictly that they obey the continuity equation ptot = ~3i/(ptot +Ptot)- 

So, one can express C as a function of (50 and therefore as a function of C, finding (after a lenghty computation) 
that C = C. 
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